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H ' Abstract 

The perturbation of the Dirac sea to first order in the external potential is calculated 
in an expansion around the light cone. It is shown that the perturbation consists of a 
causal contribution, which describes the singular behavior of the Dirac sea on the light 
cone and contains bounded line integrals over the potential and its partial derivatives, 
^ , and a noncausal contribution, which is a smooth function. As a preparatory step, we 

OO ' construct a formal solution of the inhomogeneous Klein-Gordon equation in terms of an 

04 , infinite series of line integrals. 

More generally, the method presented can be used for an explicit analysis of Feynman 
diagrams of the Dirac, Klein-Gordon, and wave equations in position space. 
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1 Introduction 

• In relativistic quantum mechanics, the problem of the unphysical negative-energy solutions 

D , of the Dirac equation is solved by the conception that all negative-energy states are occupied 

' in the vacuum forming the so-called Dirac sea. In [l], the Dirac sea was constructed for 

the Dirac equation with general interaction in terms of a formal power series in the external 
^ . potential. In the present paper, we turn our attention to a single Feynman diagram of this 

H ' perturbation expansion. More precisely, we will analyze the contribution to first order in the 

potential and derive explicit formulas for the Dirac sea in position space. Since this analysis 
does not require a detailed knowledge of the perturbation expansion for the Dirac sea, we can 
make this paper self-consistent by giving a brief introduction to the mathematical problem. 
In the vacuum, the Dirac sea is characterized by the integral over the lower mass shell 

Pi^^ y) = / + ^^p' - "^') ®(-?'°) ^""^"'"^ (1-1) 

(0 is the Heavyside function Q(x) = 1 for x > and Q(x) = otherwise); P{x,y) is a 
tempered distribution that solves the free Dirac equation (i^^ — m) P{x, y) = 0. In the case 
with interaction, the Dirac sea is accordingly described by a tempered distribution P{x, y) 
being a solution of the Dirac equation 

{i^^ + B{x)-m)P{x,y) = , (1.2) 
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where B is composed of the classical bosonic potentials. We assume ;S to be a 4 x 4 matrix 
potential satisfying the condition 7*'S(x)^7° = B{x) C^' denotes the transposed, complex 
conjugated matrix). We can thus decompose it in the form 

e = e4 + 67^ ^ + $ + i-f^E + ajk (1.3) 

with the electromagnetic potential Aj^ an axial potential Bj, scalar and pseudoscalar poten- 
tials $ and H, and a bilinear potential H^^ (see e.g. [2j for a discussion of these potentials). 
In Appendix [Bj it is shown how the results can be extended to an external gravitational field. 

The Dirac equation (jl.2p can be solved by a perturbation expansion. To first order in B, 
one has 

P{x,y) = P{x,y) + AP{x,y) + 0{B^) , 
where AP satisfies the inhomogeneous Dirac equation 

(i^^-m) AP{x,y) = -B{x) P{x,y) . (1.4) 

The factor (i^^ — m) can be inverted with a Green's function: We choose as Green's function 
the sum of the retarded and advanced Green's functions, 

six -- lun ^ f ^ + "^ e-^P(x-y) n 5) 

y> - 2 o<™ ^ J (27r)4 p^-m^± iepO ' ' ^^'^^ 

According to its definition, s satisfies the equation 

- m) s{x,y) = 6^{x-y) . (1.6) 

As a consequence, the integral 

AP{x,y) := - J d^z {s{x, z) B{z) P{z,y) + P{x, z) B{z) s{z,y)) (1.7) 
is a solution of (jl.4p . 

Clearly, p.7|) is not the only solution of the inhomogeneous Dirac equation p.4p . For 
example, we could have worked with the advanced or retarded Green's function instead of 
(|1.5p or could have omitted the second summand in p.7p . The special form of our solution 
follows from the causality principle for the Dirac sea which was introduced and discussed in 
{Xj. We will not repeat these considerations here, and simply take ()1.7p as an ad hoc formula 
for the perturbation of the Dirac sea. The reader who feels uncomfortable with this procedure 
is either referred to [T] or can, in a simplified argument, explain the special form (|1.7p from 
the "Hermiticity condition" 

AP(x,y)t = 70AP(y,x)7° , 

which seems quite natural to impose. 

In the language of Feynman diagrams, (II. 7p is a first order tree diagram. In comparison 
to diagrams of higher order or to loop diagrams, this is a very simple diagram, and it might 
seem unnecessary to study the diagram further. Unfortunately, (11. 7p gives no information on 
what AP explicitly looks like in position space. We are especially interested in the behavior 
of AP(x, y) in a neighborhood of the light cone {y — x)^ = [y — x)j{y — xy = 0. 
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Def. 1.1 A tempered distribution A{x, y) is of the order 0{{y — x)'^^), p £ Z, if the product 

{y-x)-^P A{x,y) 

is a regular distribution (i.e. a locally integrable function). It has the light-cone expansion 

oo 

Aix,y) = (1.8) 

j=9 

if the distributions A^^\x,y) are of the order 0{{y — x)"^^) and if A is approximated by the 
partial sums in the way that 

p 

^(x,y) - ^Ab](x,y) is of the order 0{{y - xfP^'^) (1.9) 

j=9 

for all p> g. 

The first summand J^^\x,y) gives the leading order of A{x,y) on the light cone. If A is 
singular on the light cone, g will be negative. Notice that the A^^^ are determined only up to 
contributions of higher order 0{{y — x)^-'"'"^), but this will not lead to any problems in the 
following. 

We point out that we do not study the convergence of the sum (jl.8p : we only make a 
statement on the approximation of A by the finite partial sums. The reason why questions of 
convergence are excluded is that the distributions A'-'^ will typically involve partial derivatives 
of order 2j of the potential B, and we can thus expect convergence only if B is analytic (for 
nonanalytic functions, the partial derivatives may increase arbitrarily fast in the order of the 
derivative, which makes convergence impossible). Analyticity of the potential, however, is 
too strong a condition for physical applications; we can assume only that B is smooth (the 
reason why analytic functions are too restrictive is that they are completely determined from 
their behavior in a small open set, which contradicts causality). Thus, the infinite sum in 
()1.8p is merely a convenient notation for the approximation by the partial sums (jl.Op . Despite 
this formal character of the sum, the light-cone expansion completely describes the behavior 
of ^(x, y) near the light cone. This situation can be seen in analogy to writing down the 
Taylor expansion for a smooth, nonanalytic function. Although the Taylor series does not 
converge in general, the Taylor polynomials give local approximations of the function. 

Our aim is to derive explicit formulas for the light-cone expansion of AP{x,y). 

2 Discussion of the Method 

Before performing the light-cone expansion, we briefly discuss the basic problem and describe 
the possible methods for calculating AP{x,y). 

At first sight, our problem seems quite complicated because of the Dirac matrices in s, 
in P, and in the potential (jl.3p . Actually, this is not the difficult point; we can reduce to a 
scalar problem by pulling all Dirac matrices out of the integral (II. 7p as follows. We have 

P{x,y) = {i^^ + m)T^2{x,y) = {-i^y + m) T^2{x,y) (2.1) 
s{x,y) = (i^^ + m) Sm2{x,y) = {-i^ + m) Sm2{x,y) , 
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where r^2 and S'„2 denote the negative-energy eigenspace and the Green's function of the 
Klein-Gordon operator, respectively: 



Tm^{x,y) = I 4^6{p'-m')Q{-p')e'^P^^~y^ (2.2) 



IS 
(27)4 

S^2{x,y) = - limyf-^^ 1 e-'P^^-y^ . (2.3) 

Using the short notation (7°)a=i,...,i6 for the bases l,i7^,7-', 7^7-', ci-''^ of the Dirac matrices, 
we can thus rewrite (|1.7p in the form 

16 

AP(x,y) = 5](% + m)7"(-i^^ + m) Ar^2[S,](:E,y) (2.4) 

a=l 

with 

Ar^2[y](x,y) 

dl^z {S^2{x,z)V{z)T^2{z,y) + T^2{x,z)V{z) Sm^iz^y)) . (2.5) 



The scalar distribution AT„2 [F](j;, y) is a solution of the inhomogeneous Klein-Gordon equa- 
tion 

(-□.-m^) Ar^2(x,y) = -V{x)T^2{x,y) , (2.6) 

as is immediately verified. Once we have derived the light-cone expansion for Ar^2(rc,y), 
the corresponding formula for lS.P{x, y) is obtained by calculating the partial derivatives and 
using the commutation rules of the Dirac matrices in (j2.4p . which will be a (lengthy but) 
straightforward computation. 

We conclude that the main problem is to calculate the solution (j2.5p of the Klein-Gordon 
equation (|2.6p . The simplest method is to analyze the partial differential equation (|2.6p . This 
hyperbolic equation is closely related to the wave equation, and the behavior near the light 
cone can be studied like the wave propagation of singularities (this method is sometimes called 
"integration along characteristics"; see e.g. [3]). In order to give an idea of the technique, we 
look at the simplified equation 

(-□ - n?) fix) = g{x) (2.7) 

and choose light-cone coordinates {u = ^(t + r), v = |(t — r), 93) around the origin (r, ■(?, Lp 
are polar coordinates in IR^). Then the D-operator takes the form 

□ = dudv - '-^{du- dv) - ^ A52 , 

where A52 = (9| + cot'd + sin~^ -i? d"^ is the spherical Laplace operator. The important 
point is that the □ -operator contains only first derivatives in both u and v. This allows us to 
express the normal derivative of / on the light cone as a line integral over / and its tangential 
derivatives. Thus, we rewrite (j2.7p on the upper light cone u = t = r, v = 0, in the form 



duiudyf (u, 0,1!), (f)) = (^du-um^ + ^As2^ f{u,0,-d,ip) - u g{u,0,'d,(p) . 
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This equation can be integrated along the hght cone as 

Jo 

= (^{du-um^ + ^As2)f{u,0,^,ip) - ug{u,0,^,(p)^ du 

= f{ui,0,^,ip) - f {0,0,^, if) 

+ J^' (^{-um^ + ^As2)f{u,0,i9,^) - ug{u,0,i3,v)J du . (2.8) 

By iterating this method, it is possible to calculate the higher derivatives in a similar way. 
We conclude that knowing / on the light cone determines all its derivatives on the light cone. 
This makes it possible to perform the light-cone expansion. We remark that complications 
arise when / has singularities on the light cone. The main disadvantage of this method is 
that the special form of the solution (12. 5p does not enter. This means, in our example, that 
additional input is needed to completely determine / on the light cone. 

Because of these problems, it is preferable to use a different method and to directly 
evaluate the integral ()2.5p . One substitutes explicit formulas for the distributions S and T 
in position space and studies the asymptotic behavior of the integral for {y — x)"^ — > 0. This 
method is presented in detail in [3]. Because it is carried out purely in position space, it gives 
a good intuition for the behavior of AP near the light cone. Unfortunately, this method is 
rather lengthy. Furthermore, the calculation of the operator products in (|2.5p and of the 
partial derivatives in (j2.4p lead to subtle analytical difficulties. 

In this paper, we use a combination of calculations in position and in momentum space, 
which gives a shorter and more systematic approach. It has the disadvantage that work- 
ing with infinite sums in momentum space is more abstract than studying the behavior of 
distributions in position space. Therefore the reader may find it instructive to compare the 
technique of this paper with the calculations in [3] . 



3 The Formal Light- Cone Expansion of ATj„2 

In this section, we will perform the light-cone expansion for AT^2{x,y) on a formal level. 
The analytic justification for the expansion is postponed to the next section. We assume that 
m / and set a = m?. 

Since we want to derive formulas in position space, it is useful first to consider explicitly 
how T^2 looks like. Calculating the Fourier transform of the lower mass shell (12. 2 p yields an 
expression containing the Bessel functions Ji, Yi, and Ki. The most convenient form for our 
purpose is to work with the power series for these Bessel functions, which gives 

Ta{x,y) = --^ lim 



87r3 0<e^0 ^2 _ ^^^0 

327r3 o<e^oV + 



EinTi^^r^m + i) + m) (3.1) 



327r3 ^^11(1 + 1)! 4' 
with ^ = {y — x), c = 2C — log 2 with Euler's constant C, and the function 

n 

$(0) = , $(n) = Y^T for n > 1 . 

k=i ^ 
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The logarithm is understood in the complex plane, which is cut along the positive real axis 
(so that limo<e^o log(a; + ie) = log|x| is real for x > 0). It can be verified explicitly that 
Ta is a solution of the Klein-Gordon equation {—Dx — a) Ta{x,y) = 0. Furthermore, one can 
calculate the Fourier transform Ta{p) with contour integrals. For > 0, the ^*^-integral can 
be closed in the lower complex plane, which gives zero. In this way, one immediately verifies 
that Ta is formed only of negative-energy states. 

This formula for Ta looks quite complicated, and we do not need the details in this section. 
It suffices to observe that Ta has singularities on the light cone of the form of a pole and a 
^-distribution, 

1 PP 

where PP denotes the principal value. Furthermore, there are logarithmic and 0-like contri- 
butions, since 

lim log{af-ief) + in = log{\aC^\) + in Q{^^) e{f) , 

0<e^0 

where e is the step function e(x) = 1 for x > and e{x) = — 1 otherwise. The important 
point for the following is the qualitative observation that the contributions of higher order 
in a contain more factors and are thus of higher order on the light cone. This yields 
the possibility of performing the light-cone expansion by expressing ATa in terms of the 
a-derivatives of Ta- In the following lemma, we combine this idea with the fact that line 
integrals over the potential should occur according to (j2.8p . The lemma gives an explicit 
solution of the inhomogeneous Klein-Gordon equation (j2.6p and is the key for the light-cone 
expansion. We use the notation 



= (-] T 



da 



a 



Lemma 3.1 The formal series 

A{x,y) = («-«'r(n"V^)|a,+(i-a)xdari"+i)(x,y) (3.2) 

n=0 "'O 

satisfies the equation 

(-□, - a) A{x, y) = -V{x) Ta{x, y) . (3.3) 
Proof: In momentum space, Ta has the form 

Ta{p) = Hp^ - a) e(V) • 

Since a > 0, the mass shell does not intersect the hyperplane p^ = 0, and we can thus 
calculate the distributional derivative to 

^ (jW (p2 _ a) e(-p°)) = 2pj (^("+1) (p2 _ a) e(-p°) - (p2 - a) <5° 5(p°) 

= 2pjj("+i)(/_a)e(V) • 
Hence, for the calculation of derivatives we can view Ta as a function of (p^ — a); i.e. 

At(")(p) = -2p,r("+i)(p) . (3.4) 
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(n) 

This relation can also be used to calculate the derivatives of Ta in position space, 



7 (27r)4 



(2vr)^ 

2 y (27r)4 5pi 

2 y (27r)4 ^pj "^ 

= ^(y-x),Ti"-i)(x,y) . (3.5) 

Using that Tq is a solution of the Klein-Gordon equation, we also have 

= y-£) {p'-a)Up) = {p'-a)T^-\p) - nTj,--'Hp) 

and thus 

(-□.-a)rW(x,y) = nr("-i)(x,y) . (3.6) 



With the help of (jS.Sp and (j3.6p . we can calculate the derivatives of the individual summands 
in (13.21) as follows: 



(-□. - a) [\a - a^r (□"^)|„j,+(i-a). da ri"+i)(x,y) 

\l _ «)2 (a - (□"+V)|,,+(i_,), da ri"+i)(x,y) 
\l -a) (a- a^r (5,n"l^)|«j,+(i-a). da {y - ri")(x,y) 
+(n + l) /'(a-a2)"(n"y)|,,+(i_„),darW(x,y) . 



In the second summand, we rewrite the partial derivative as a derivative with respect to a, 

and integrate by parts, which gives 

(-□x - a) [\a - a^T (□"V^)|ay+(i-a)x da Ti"+'\x,y) 

•J 

= 5nfiV{x)Ta{x,y) 

+nj\l - af (a - a^)""! (□"y)|,,+(i_,), da Ti^\x,y) 

- f\l - af (a - a^r (□"+V)|,,+(i_,), da T^^^'\x,y) . 

^ 

After dividing by n! and summation over n, the last two summands are telescopic and vanish. 
This yields ■ 



It would be nice if the solution of the inhomogeneous Klein-Gordon equation constructed 
in the previous lemma coincided with ATa{x,y). This is really the case, although it is not 
obvious. In the remainder of this section, we will prove it. The technique is to expand (j2.5p 
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in momentum space and to show that the resulting expression is the Fourier transform of 
([32]). 

Since (|2.5p is Hnear in V, we can assume that V has the form of a plane wave, 

V{x) = e"^"^ . (3.7) 
Transforming ATa{x,y) to momentum space gives the formula 
J (fx J (fy ATa{x,y) e^^^jz-ipix 

PP \ 

Ta{pi) + Ta{p2) 6\q-p2+pi) , (3.8) 

pi - a J 



PP 

pl-a ■ ~ ' Pt 

where pi,P2 are the in- and outgoing momenta, and where PP denotes the principal value, 

± 



— — hm 

p-^ — a 2 o<£^o ^ p^ — a± lep^ 



The factor 5^{q — P2+P1) in (j3.8p describes the conservation of energy- momentum. In order 
to simplify the notation, we leave out this 6^-factoT and view ATq as a function of only one 
free variable p = {pi + P2)/2, 

q q PP q g PP 

AUp+-,p--) := - (^p^^y_J aiP-^)-Ta{p+^) _ _ ^ ■ (3-9) 



The transformation to position space is then given by a single integral, 

Ar„(x,y) = y'^Ar,(p + |,p-|)e-f(--^')e-i(--+^) , (3.10) 
the inverse transformation to momentum space takes the form 

ATa{p + \,p-\) = Jd'yATa{x,y)e'P'^^-y^e^i^-+y^ . (3.11) 

The first step for the light-cone expansion in momentum space is to rewrite ()3.9p in the 
form 

^Ta{p+ -,P- -) = -77 ^ 1 , ^ Ta{p--) - Ta{p + 



2'^ 2' [(p_|)2_a] + 2pg -^^^ 2' '^^^ ' 2' [{p + if -a]-2pq 

and to use that the expressions within brackets [• • •] vanish as the arguments of the 5- 
distributions in Ta(p it |), 

AT.(p+|,p-|) = g(r.(p+|)-T.(p-|: 

Now we expand Tq in a Taylor series in q, 

AUp + -,P--) = — 2 g ^^^^ ^^^^^^^ r.(p) - • • • ^ . (3.12) 

We want to rewrite the ^-derivatives as derivatives with respect to a. This can be done by 
iteratively carrying out the p-derivatives with the differentiation rule (j3.4p . One must keep 
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in mind that the p-derivatives act either on or on the factors pj that were previously 
generated by (j3.4p . e.g., 



^2 



In this way, carrying out the p-derivatives in (j3.12p gives a sum of many terms. The combina- 
torics may be described as follows: Each factor that is generated by ()3.4p and differentiated 
thereafter gives a pairing between two of the derivatives , . . . , di^^.^^ ; namely, between the 
derivative by which it was created and the derivative by which it was subsequently annihi- 
lated. The individual expressions obtained after carrying out all the derivatives correspond to 
the possible configurations of the pairings among the di^ - ■ ■ di^^.^^ ■ They depend only on the 
number of pairs and not on their specific configuration. More precisely, every pair increases 
the degree of the derivative of Ta^ by one and gives a factor q'^/2, whereas the unpaired 
derivatives also increase the degree of Ta^ and give a factor pq. 

It remains to count how many configurations of such n pairs exist. We use the notation 

m 
n 

combinatorics becomes clearer if one first selects 2n out of the m points and then counts the 
number of possible pairings among these 2n points to (2n — 1)!!. This explains the formula 



for the number of possibilities to choose n pairs from a set of m > 2n points. The 



m 




( m\ 


n 




\2n) 



TTl ' 1 

(2n-l)!! = -— . (3.13) 

^ ^ (m - 2n)! 2" n! ^ ^ 



We conclude that 



Q2k+1 qii qi2k+l 

-Ta{p) 



dph . . . dp'2k+i "-^^'2 2 

= Vf (2fc + l)! 1 ^(2fc+l-n)/ X / x2fe+l-2n 

' {2k + l-2n)\2^n\ " V^/ 
After substituting into (j3.12p and reordering the sums, we obtain 

Finally, we pull one factor pq out of the sum, which cancels the principal value, 

AT„(.. |..- 1) ^ -|tir (^)"|_M!^^..««,(„ ^ (,15, 

This is the formula for the light-cone expansion in momentum space. 

It remains to show that the Fourier transform (j3.10p of (j3.15p coincides with (|3.2p . In 
order to see a first similarity between these formulas, we substitute the plane-wave ansatz 
(13771) into 



71=0 " 2 



(3.16) 



The last exponential factor also occurs in (|3.10p . Furthermore, it is encouraging that both 
()3.15p and ()3.16p contain a power series in q^. The main difference between the formulas 
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is related to the factor exp(—iTq{y — x)) in (j3.16p : Expanding this exponential yields a 
power series in q{y — x). In momentum space, this corresponds to a power series in q^dpj 
(because differentiating the exponential in (|3.10p with respect to p gives a factor i{y — x)). 
The expansion (|3.15p . however, contains a power series in pq, not in q-'dpj. The following 
lemma allows us to transform these expansions into each other. 



Lemma 3.2 For all r > 0, 

2k 



q^ d 



Ti'Hp) = E(-i)' 



1=0 
k 



2k 
I 



{pqr T(^^''Kp) = E 
1=0 



2k 
I 



qi d 
'2dp^ 



(pg)2fe-2«^^(r+2fc-0(p) (3.17) 
2k-2l 



(3.18) 



Proof: On the left side of (j3.17p . we calculate the derivatives inductively using the dif- 
ferentiation rule (|3.4p . The derivatives either act on Ta \ which increases the order of the 
derivative of Ta^ and generates a factor pq, or they act on previously generated factors pq, 
which reduces the number of factors pq by one and produces a factor q'^/2. This can also be 
written in the inductive form 



-pq 



q^ d 



l~2 



2 dpi 



2 dp) 



The combinatorics is described by counting the number of possibilities in forming / pairs 
among the 2k derivatives. 

Equation ()3.18p follows in the same way from the relation 



2 dp^ 



(P) ■ 



After these preparations, we can prove the main result of this section: 



Theorem 3.3 (formal light-cone expansion of AT„j2) For m ^ 0, the distribution 

AT„2(x,y), i2. 5\) . has a representation as the formal series 



Ar^2(x,y) = («-«')"(n"^)|ay+(i-a).daT^"+')(x,y) 

n=0 "'O 



(3.19) 



Proof: We expand the factor (^ — r^)" and the exponential exp(— irg(?/ — x)) in (j3.16p . and 
so have 



A{x,y) 



n=0 



.2^n 4 n 



X 

2 l=Q 



n 



^2\l 



1\"-' 

4. 



^„ n! A: 



n=0 



X e-'^'^^y-''^ dr ri"+^)(2;, y) e'^a i^+v) 

n 



k=0 1=0 



X / r'=+2^dTri"+i)(x,y)e-*§(^ 



+y) 



(3.20) 
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Next we carry out the r-integration. This gives a contribution only for k even, 



^ i-^ ^ i-iq{y - x)r V, . y 
^ ^ {2k)\ ' 



n=0 

oo 

E- 

n=0 



n! 



fc=0 



/=0 



n! 



A i ^ # (2A;)! ^ 



1 

4 



n—l 



k=0 



1=0 



2k + 21 + I 



{2k + 21 + 1) 4*=+' 

ri"+i)(x,y)e-*5(-+2/) 



Now we transform to momentum space by substituting into (j3.1ip . The factors —iq{y — x)/2 
can be rewritten as derivatives ^d^j acting on the plane wave e^^^^~y\ Integrating these 
p-derivatives by parts gives 



Aip+l^p-^) 



E 

?i=0 



(-1)" 



fc=o 



'-Xf—y. 



-ly 



2k 



n 



1=0 



2k + 21 + 1 \l 



We shift the indices n and k according to n — / — > n and k + I 
of the ^-summation to / = 0, . . . , A;. We thus obtain 



E- 

n=0 

oo 

E- 

n=0 



n 



k=0 

oo 



2fc + 1 (2A; - 2/)! /! \^ 4 



n 



A:=0 



{2k + 



1)1 ^ 



2k 
I 



k. This changes the range 



2fc-2/ 



2 Sp- 



in this way, we have transformed the line integrals in (j3.16p to momentum space. We remark 
that we did not use the special form of Ta in the calculation from (I3.20p to (13.210 : for the 
computation so far, we could replace Ta by any other function. 

In the last step, we apply Lemma [32} substituting ()3.18p into the light-cone expansion 
(j3.15p for ATa also yields the expression (|3.2ip . Thus A coincides with ATq, which concludes 
the proof. ■ 



4 Resummation of the Noncausal Contribution 

In this section, we will put the previous formal calculations on a rigorous basis. The in- 
teresting part is to recover the noncausal structure of AT^2(3;,y) by resumming the formal 
light-cone expansion. We begin with specifying the conditions on the potential V in ()2.5p . 

Lemma 4.1 Let V G L^(IR'^) be a potential which decays so fast at infinity that the functions 
X* V{x) are also . Then AT^2{x,y) of 112. 5\) is a well-defined tempered distribution on 
IR4 X IRl 

Proof: It is easier to proceed in momentum space and to show that 

^T^2{p2,Pl) = -S^2{p2) V{p2 -Pl) T^2{pi) - r„2(p2) V{p2-Pl) Sm2{Pl) (4.1) 

is a well-defined tempered distribution, where V is the Fourier transform of V. The assump- 
tion then follows by Fourier transformation. 
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In momentum space, the conditions on the potential give V £ C1(IR4) n L°°(IR'^). We 
choose two test functions f,g£ 5(IR^). Then the function V{p2 — •) /(•) is and has rapid 
decay at infinity. Thus the integral over the lower mass shell 

is finite and depends differentiably on p2. Consequently, the product gl is in and has 
rapid decay, and we can calculate the principal value by 



^(^^2) := / j—rgV{p2-pi)T^2{pi)f{pi) 



which gives a finite number. This shows that the first summand 

Sm-i {P2) V{p2 - Pi) T„2 (pi) 

in (j4.ip is a well-defined linear functional on 5'(IR'^) x S'(IR^). This functional is bounded in 
the Schwartz norms ||.||o,0) ll-lk.O; IMIo,i) IMU,! of f,g, which gives continuity. 

For the second summand in (j4.ip , one can argue in the same way after exchanging pi and 
P2- ■ 

For the light-cone expansion, we clearly need a smooth potential. Therefore, we will assume 
in the following that V £ C°° n and V{x) £ L^. 

We come to the mathematical analysis of the light-cone expansion. We again assume 
that m 7^ 0; the case m = will be obtained at the end of this section in the limit m ^ 0. 
In the first step, we disregard the convergence of the infinite sums and check that all the 
performed operations make sense and that all expressions are well-defined. We start with 
the end formula (13.19P of the light-cone expansion. The line integrals over the potentials 
are C°*^-functions in x,y. The factors T^}{x,y) are tempered distributions, as one sees after 
differentiating the explicit formula (j3.ip with respect to a. Thus (j3.19p makes mathematical 
sense. The calculations leading to this result are not problematic except for the handling of 
the principal value following (j3.9p . The easiest method for studying this more rigorously is 
to regularize the principal value in (j3.9p with the replacement 



PP 1 1 



X 2 X zizie + 

Then all the subsequent transformations are well-defined, and the critical operation is the 
cancellation of the principal value against one factor pq before (j3.15p . In order to justify this 
operation, we use in (|3.14p the exact formula 



PL_ fpqfk+i ^ {iPQ? + e')-e' 2k ^ £! | . .2k 

{pqY + £2 ^^^^ (pg)2 + ^2 ^^P^) y (p^)2 + £2 I ^^P^) 



2k+2 

{pqf^ - (pqf^-^ + ••• + (-1)'= e^fc + ^ 



{pqy + 

The first summand gives (j3.15p : the following summands (—1)' e^' {pq)'^'^~'^^ contain no prin- 
cipal value and vanish in the limit e — > 0. Thus it remains to consider the last summand for 
e ^ 0, 



g2fc+2 



{pqY + 



^(2.-M+.)(^) . (4.2) 
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We use that the support of T^2{p) is on the mass sheh = m? and apply the relation 



£^o(pq')2 + e2 a. \daj e^Q {pqY + e"^ 



da / 

This expression is well-defined for m 7^ 0. The limit (j4.2p contains an additional factor 
e^'^+^((jf^/4)" and thus vanishes. 

We conclude that the light-cone expansion is mathematically rigorous except for the 
formal character of the infinite sums. In the remainder of this section, we will carefully 
analyze the infinite sum in (j3.19p . More precisely, we will do the following: As explained in 
the introduction, the infinite sum in (jl.Sp is only a notation for the approximation (11.90 of the 
partial sums. Following this definition, we need only show that the light-cone expansion is 
well-defined to any order on the light cone, but we need not study the convergence of the sum 
over the order on the light cone. According to the explicit formula (13. ip . each factor T^"*"^^ 
in (j3.19p consists of an infinite number of terms of different order on the light cone (we will 
see this in more detail in a moment). In order to bring (j3.19p into the required form (jl.Sp . 
we must collect all contributions to a given order on the light cone and form their sum. This 
procedure is called resummation of the light-cone expansion. If the sum over all contributions 
to a given order on the light cone were finite, this resummation would be trivial; it would 
just correspond to a rearrangement of the summands. It will turn out, however, that these 
sums are infinite, and we must find a way to carry them out. 

In order to see the basic problem in more detail, we consider the explicit formula (|3.ip 
for Ta- We start with the last summand 

n °° (--1)1 (np'^V 
^■.(-.rf - -S^i^g^ 4^ + (4.3) 

(the notation 'x' means that we consider only a certain contribution to Ta). This is a power 
series in a, and we can calculate its derivatives to 

- i^E 1 w + + . »>i. (4.4) 

l=n—l ^ ' 

For increasing n, the derivatives are of higher order on the light cone; more precisely, the 
contribution (j4.4p is of the order 0((y — x)^^"^^)). Thus, the contribution of (|4.3p to the 
formal light-cone expansion ()3.19p consists, to any order on the light cone, of only a finite 
number of terms. Thus the resummation is trivial. Of course, we could rearrange the sum 
by collecting all the summands in (|4.3p and (j4.4p of a given degree in ^ and writing them 
together, but this is only a matter of taste and is not really needed. 
For the second summand in (13.11). 



U..y) X -^(M{^-.iO«») + i, + .)^^i^M_i (4.5) 



1=0 
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the situation is more complicated: The contribution (|4.5p is a power series in a and can 
be discussed exactly as ()4.3p . The only difference (apart from the missing factors $) is the 
prefactor (log(^^ — iOS,^) + zvr + c), which has a logarithmic pole on the light cone. The 
contribution ()4.6p . however, contains a factor logo and is not a power series in a. As a 
consequence, the higher a-derivatives of (j4.6p are not of higher order on the light cone. For 
example, the contribution to the order 0{{y — x)'^) has the form 

Ta{x,y) X a log a + 0{{y - x f) 

Ti^Hx.y) X -^{l + \oga) + 0{{y-xf) 

(_-[\n 1 

rW(x,y) X L_L(„-2)!-^ + , n>2 . 

This means that we must resum an infinite number of terms; more precisely, 

1 /"^ 

AT]0](x,y) X y\ay+{i-c.)x da {I + logo) 

+1^ E - / (« - ^« + ■■■ + 0{{y - xf) . (4.7) 

n=l ^ 

This is a serious problem. Namely, we can expect the series in ()4.7p to converge only if the 
derivatives d^V do not grow too fast in the order 2n of the derivative. It turns out that 
analyticity of V is necessary for convergence, which is too restrictive. 

On a technical level, this convergence problem of the contributions to Ar^2 to a given 
order on the light cone is a consequence of the factor log a in (j4.6p ; we call it the logarith- 
mic mass problem. Because AT^2(x,y) is well-defined by (12. 4p . it is not a problem of the 
perturbation expansion, but rather shows that the light-cone expansion was not performed 
properly. The deeper reason for the convergence problem is that we expressed Ar^2(x,y) 
only in terms of the potential and its derivatives along the line segment xy. However, the 
perturbation AT^2{x,y) is not causal in this sense; it depends on V in the whole Minkowski 
space (this becomes clear in ()2.4p from the fact that the support of T^2{x,.) is IR^). In 
a formal expansion, we can express AT^2(x,y) in terms of □'^l^Aj/+(i-A)a;) < A < 1, but 
we cannot expect this expansion to converge. The simplest one-dimensional analog of this 
situation is the formal Taylor series 

oo 1 

1 1 . 

of a smooth function. The right side cannot in general converge, because it is not possible to 
express /(x), x 7^ 0, in terms of /^"^(O). 

The solution to the logarithmic mass problem is to reformulate the problematic contribu- 
tion of ()4.6p to the light-cone expansion (13.19P as a noncausal term that is obviously finite. In 
some sense, we will simply reverse our former construction of the light-cone expansion. This 
is not trivial, however, because the differentiation rule (13. 4p . which was crucial for rewriting 
the Taylor expansion (j3.12p as an expansion in the mass parameter a, is not valid for (j4.6p . 
In the end, we want to write the light-cone expansion in a way which shows that part of the 
behavior of AT„2(x,y) can be described with line integrals of the form (I3.19P whereas other 
contributions are noncausal in a specific way. 
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We work in momentum space. The Fourier transform of the problematic series (|4.6p is 

Up) = /.*....^_LgJ^^ = . (4.8) 

Notice that this expression is highly singular at p = 0; especially, it is not a distribution. 
However, it is well-defined as a distribution on analytic functions in p. This comprises all 
functions with compact support in position space, which is a sufficiently large function space 
for the following. Furthermore, we introduce the series 

oo I 

Laip) = ^E^o'S'ip) (4.9) 



and set 

j(n) _ /^A^" - '^"^ / d ^ " 

'^'^ " [da. 



Lemma 4.2 The series ^4-^ o,nd satisfy the relations 

Ja{p) = [\ra{p)dT (4.10) 



^L(^Hp) = -2pjL(^^'Hp) . (4.11) 

Proof: Equation ()4.10p is verified by integrating the power series (14. 9p and comparing with 
()4.8p . The distribution pj 5'^{p) vanishes identically. Since the derivatives of distributions are 
defined in the weak sense, it follows that 

= □"+i(p,<54(p)) = p,n^^'6\p) + 2in + l)^n^6\p) , 

and thus 

Applying this relation to every term of the series dM]) yields that dp, La = -2pjL^a \ and 
(j4.1ip follows by differentiating with respect to a. ■ 

The function La is useful because (14. lip coincides with the differentiation rule (|3.4p for T^. 
This implies that all the formulas for Tq, especially the manipulations of Lemma [3.2[ are also 
valid for La- 

The following technical lemma is the key for handling the logarithmic mass problem. 
Lemma 4.3 (resummation of the noncausal contribution) // V is the plane wave 



oo , „i / d\ 

- E (« - da (^— j (a log(a) Ja{p)) (4.12) 

[ dfi ( {a + b) log{a + b) [ L^a+iT~i)b~f,pq{p) dr] . (4.13) 



n=0 

_1 d_ 

~2 da 
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Proof: The series (|4.12|) is obtained from the formula (j3.2|) for A{x,y), by the replacement 
Ta — > alog{a)Ja. As remarked in the proof of Theorem 13.31 ^-H t^is transformations from 
(|3.20p to (|3.2ip are also valid if we replace by any other function. Therefore, 



Km 



n=0 



nl 



[ 2 



2k-2l 



2k 
I 



d 
da 



n+l+l 



(a log(a) Ja{p)) ■ 



We carry out the sum over n by redefining a as a = m —q /4 and substitute (j4.10p as follows: 



2A; 
/ 



d 

'J dpi 



2k -21 



da 



j log (a) ^ Lraip) dr^ . 



Using that La{p) and Ta{p) obey the analogous differentiation rules (|4.1ip and (j3.4p . respec- 
tively, we can apply relation (j3.17p with Ta''^ replaced by lI'^'' to obtain 



2A; 



fe-/ 



s=0 



2A; - 2/ 
s 



l+s 



x{pq) 



2k-2l-2s 



da 



i+i 



log(a) / L^^t^'-^\p)dT 



We introduce the index r = 1 + s, replace s by r — / , and substitute the combinatorial formula 
(|3T3]) : 



(liT2D 



2fc + 1 ^ r! 



1 



A;=0 



r=0 

r 

xE(-i) 

/=0 



{2k - 2r)\ \ 4 ^ 



2k-2r 



I I r 

I I \ da 



The last sum can be eliminated using the combinatorics of the product rule, 



E(-i 

1=0 



* W ) (da 



(a log(a)4i'-2^^-+(-'))(p)) 

Furthermore, we shift the index k according to k — r ^ k, yielding 

1 



KT2\i 



E 



2-^ ^1 /I 



j42t + 2r + l (2*:)! r! \^ 4 ^ 
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Without the factor {2k + 2r + l)^"*^, we had two separate Taylor series which could easily be 
carried out explicitly. The coupling of the two series by this factor can be described with an 
additional line integral, 



2r 



X 



d_ 

jib 

We finally carry out the remaining Taylor sum 



a + b) log(a + b) f L^a+{r-i)b+miP) '^'^ 

•J 



|b=0 



The result of this lemma is quite complicated. The important point is that the convergence 
problems of the infinite series ()4.12p have disappeared in (14.13p . which is obviously finite. 
Namely, the a-derivative of the integrand in (j4.13p has at most logarithmic singularities. 
These singularities are integrable and disappear when the //-integration is carried out. 

After these preparations, we can state the main theorem. Since the resulting expansion 
is regular in the limit m ^ 0, it is also valid for m = 0. 

Theorem 4.4 (light-cone expansion of Ar^2) The distribution Ar^2 of 12. 5\) has the 

representation 

AT^2{x,y) = / (a-ay (□"y)|„,+(i-.).darj("+^)(x,y) (4.14) 

+N^2{x,y) (4.15) 



with 

TP = Ta- a log(a) J, , T^ in) ^ (A. ^ ^reg 



da, 

and the Bessel function Ja of |^.<§[ ). The series j^[ ) is well-defined in the sense of Def. 

The contribution Nj^2{x,y) is a smooth function in x,y and has a representation as the 
Fourier integral 

N^2{x,y) = J^J^V{q)N^2ip,q)e-^P^^-y^e-^i^^+y^ wzth (4.16) 

Na{p,q) = -\^f d, log fa - (1 - ^) /3 Mp)f--it';' . (4.17) 
I daj-i y 4y i3=-fi^^+fipq 

Proof: By definition, T^^^ differs from Ta by the contribution (j4.6p . Thus an explicit formula 
for T^^^ is obtained from (13. ip if we replace the factor log(a^^ — ie^^) in the second line by 
log(^^— ie^^). As a consequence, T^^^ is a power series in a, and the higher order contributions 
in a are of higher order on the light cone. This justifies the infinite sum in (|4.14[) in the sense 
of Def. 11.11 Furthermore, the difference between the formal light-cone expansions (j3.19p and 
(j4.14p coincides with the contribution (I4.12p . which was resummed in Lemma 14.31 We carry 
out the r-integration in (j4.13p using the series expansions (j4.8p and (j4.9p , which gives (j4.17p . 

The g-integral in (I4.16P is well-defined since V{q) is and decays sufficiently fast at 
infinity. Finally, the p-integration can be carried out with the ^^-distributions in ()4.8p , which 
gives a smooth function N^2{x,y). H 

We call (|4.14p and (j4.15p the causal and noncausal contributions, respectively. 

We could proceed by studying the noncausal contribution more explicitly in position 
space. For the purpose of this paper, however, it is sufficient to notice that N{x, y) is smooth 
on the light cone. 
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5 The Light-Cone Expansion of the Dirac Sea 



Having performed the light-cone expansion for Ar^2, we now return to the study of the 
Dirac sea (jl.7p . From the theoretical point of view, the light-cone expansion for AP^2 is 
an immediate consequence of Theorem 14.41 and formula (|2.4p : We substitute the light-cone 
expansion (|4.14P " (j4.15p into (|2.4p . Calculating the partial derivatives and dy of the causal 
contribution (j4.14p gives expressions of the form 

AP^2{x,y) X /V(«)D»n%,+(i_A).daOTj("+^)(x,y) , (5.1) 



which are again causal in the sense that they depend on the potential and its partial deriva- 
tives only along the line segment xy (here V{a) denotes a polynomial in a; D"^ stands for any 
partial derivatives of the order a). Since (|5.ip contains distributional derivatives of T^f^"'~^^\ 
it is in general more singular on the light cone than the corresponding contribution to Ar^2 . 
On the other hand, the partial derivatives of the noncausal contribution N^2{x,y) can be 
calculated with (14.171) and yield smooth functions in x, y. We conclude that the qualitative 
picture of Theorem 14.41 especially the splitting into a causal and a noncausal contributions, 
is also valid for the Dirac sea. 

The situation becomes more complicated if one wants to go beyond this qualitative picture 
and is interested in explicit formulas for the Dirac sea. The problem is to find an effective 
and reliable method for calculating the partial derivatives and for handling the combinatorics 
of the Dirac matrices. Before entering these computational details, we explain how the 
qualitative picture of the light-cone expansion can be understood directly from the integral 
formula (jl.Tp : The tempered distributions s{x,y) and P{x,y) are regular functions for (y — 
x)'^ 7^ and are singular on the light cone (this can be seen explicitly from e.g. (12. ip and 
(|3.ip l. Integrals of the form 

J P{x, z) f{z) (fz or j s{x, z) f{z) (fz 

with a smooth function / (which decays sufficiently fast at infinity) give smooth functions 
in X. The integral in (jl.Tp is more complicated because it contains two distributional factors 
s,P. This causes complications only if the singularities of s and P meet; i.e., if z lies on the 
intersection Lx H Ly of the light cones around x, y, where 

= {y G IR\ (y - xf = o} . 

If y—x is timelike or spacelike, then LxCiLy is a 2-sphere or a hyperboloid (respectively), either 
of which depends smoothly on x, y. As a consequence, the integral over these singularities can 
be carried out in ()1.7p and gives a smooth function. On the light cone (y — x)^ = 0, however, 
Lx n Ly does not depend smoothly on x, y. More precisely, in the limit < (y — x)^ — > 0, 
the 2-sphere L.^ n Ly degenerates to the line segment {Ay + (1 — Ax) , < A < 1}. The limit 
> (y — x)^ 0, on the other hand, gives the degenerated hyperboloid {Ay + (1 — Ax) , A < 
or A > 1}. This simple consideration explains why the singularities of AP(x,y) occur on 
the light cone and makes it plausible that the behavior of the singularities is characterized 
by the potential and its derivatives along the line xy = {Ay + (1 — A)x , A G IR}. Clearly, 
V{z) also enters into AP^2(x,y) for z xy, but this noncausal contribution is not related to 
the discontinuity of L^ Ci Ly on the light cone and is therefore smooth. The special form of 
the singularities, 

AP(x,y) ~ D"log((y-x)2-i0(y-x)°) (y-x)2" , (5.2) 
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is less obvious. That the potential enters only along the line segment xy can be understood 
only from the special form of (jl.7p : it is a consequence of the causality principle for the Dirac 
sea that was introduced in [1]. In fact, it gives an easy way to understand the meaning of 
"causality" of the perturbation expansion for the Dirac sea. 

We finally describe our method for explicitly calculating AP{x,y). As in Theorem 14.41 
we will not study the noncausal contribution; we are content with the fact that it is bounded 
and smooth. In other words, we consider only the singular contribution (jS.ip to the Dirac 
sea. Since the difference between T^"'} and T^^^ is smooth, we can iust as well consider 
the formal light-cone expansion (j3.19p and calculate modulo smooth terms on the light cone. 
This has the advantage that we can work with the useful differentiation rule ()3.5p . The 
calculation can be split into several steps, which may be listed as follows. 

(1) Calculation of the partial derivatives with the product rule and the differentiation for- 
mulas 



d 
dyi 



(5.3) 

^V(a)i?'^n = ^'aP(a)9,Z)«n%,+(i_„), (5.4) 

^ J^' Via) D^u'V^^y+ii.^), = ^V-«)^(«)fi'j^"n'W{i-a)x (5.5) 
^(y-2;)fc = —^{y-x)k = -Qjk ■ (5.6) 

(2) Simplification of the Dirac matrices with the anti-commutation relations {7-^ , 7^^'} = 2g^^ . 
This leads to a contraction of tensor indices. The generated factors {y — x)^ and 
(y — xydjV are simplified in the following calculation steps (3) and (4). 



(3) Absorption of the factors {y — x)^. We calculate the Laplacian by iterating (|5.3p . 

□ .Tir^)(x,y) = -2Tir^)(x,y) + i(y-x)2ri'l)(x,y) , 
and then combine it with (j3.6p . which gives the rule 

{y-xfT^^Hx,y) = -4nr^"+^)(x,y) - Am^ T^^^^\x , y) . (5.7) 

(4) Partial integration of the tangential derivatives^ 

^V(a)(y-x)^9,D"n%j,+(i_,), = j\{a) -^D-n'V\^y+^^_^), 



P(a)Z)'^n%,+(i_,w - / 7''(a)I?'^n . (5.8) 



0=0 Jq 

After these steps, AP^2(x,y) consists of many terms of the form 

AP^2(x,y) X (causal expression in D°-\2^V ) x T^^}{x,y) , n > — 1 
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It remains to insert the series representations for T^"^(x,?/). It is useful first to introduce the 
short notation = c^^" (n > 0) and 



Z-2 



— hm — 

2 0<e-^0^^2_^g^0 



0\2 



log 2: := ^ ^^hn ^ ^^Sif - 
(5) Substitution of the explicit formulas 



(-1). ^ I o . 1 m2'+4 , ^ ^ 

{x,y) X " ^ ^^^ir(7T2)!^ ^^-^^ 



I 


- 


2 

m 


1 




8^ 


1 




1 






8^ 


1 


00 

E - 

Z=n-1 




+1 ^2i+2-i 




4^+1 l\{l + l- 



Z=n— 1 ^ ' 

where again we have used the notation of (13. ip (we take only the singular contribution 
on the light cone; T^-^'^ is defined via (|5.3p ). 

In this way, the calculation of the causal contribution is reduced to a small number of symbolic 
computation rules (I5.3p - (l5.1ip . which can be applied mechanically. This makes it possible 
to use a computer program for the calculation. The C++ program "class_commute" was 
designed for this task (commented source code available from the author on request). It 
computes the causal contribution for a general perturbation (II. 3p to any order on the light 
cone. The formulas to the order 0{{y — x)^) modulo the noncausal contribution are listed in 
the appendix. 



6 Outlook 

In this paper, the light-cone expansion was performed for the Dirac sea to first order in the 
external potential. The presented method can be generalized in several directions and applied 
to related problems, which we now briefly outline. 

First of all, the method is not restricted to the Dirac and Klein-Gordon equations; it can 
also be used for the analysis of other scalar and matrix hyperbolic equations (in any space- 
time dimension). The consideration (|2.8p . which gives the basic explanation for the line 
integrals in the light-cone expansion, can be applied to any hyperbolic equation (in curved 
space-time, the line integrals must be replaced by integrals along null geodesies; see e.g. [3]). 
Thus, the behavior of the solution near the light cone is again described by an infinite series 
of line integrals. The line integrals might be unbounded, however, which leads to additional 
convergence problems (e.g., one can replace the integrals in (j3.2p by ^ /f^ e(a) do • • •, which 
gives a different formal solution of (j3.3p ). 

Furthermore, the light-cone expansion can be generalized to symmetric eigensolutions and 
the fundamental solutions: The formal light-cone expansion of Section 3 applies in the same 
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way to any Lorentzian invariant family T„2 of solutions of the Klein-Gordon equation, i.e. to 
a linear combination of 

T^2{p) = Sip^-m^) and (6.12) 

Tm<p) = Sip" - m") eijP) . (6.13) 

The second case ()6.13p allows us to generalize the light-cone expansion to the Green's function. 
The advanced Green's function of the Klein-Gordon operator, for example, can be derived 
from T„2 of (|6.13p by 

5^2 (x,y) = 27rir^2(x,y)G(yO-xO) . 

This relation even remains valid in the perturbation expansion, e.g. to first order, 

/^Sl,2{x,y) = 2^iAT^2(x,y)e(y°-x°) (6.14) 

(for a derivation of this formula in the context of the Dirac equation, see [1]). Thus the 
light-cone expansion for AT^2 immediately yields corresponding formulas for the Green's 
function. 

In contrast to the formal light-cone expansion of Section 2, the resummation of the non- 
causal contribution depends much on the particular problem. An analysis in position space 
according to [4J might be helpful for the understanding of the noncausality. For X^2 accord- 
ing to (I6.13p . for example, there is no noncausal contribution at all, which also simplifies the 
analysis of the Green's functions. 

By iteration, the method can also be applied to higher-order Feynman diagrams and even 
makes it possible to sum up certain classes of Feynman diagrams explicitly. For the Dirac 
Green's function and the Dirac sea, this is explained in detail in [5]. 

A Some Formulas of the Light-Cone Expansion 

The following formulas give AP(x, y) to first order in the external potential (|1.3p up to 
contributions of the order 0{{y — x)^) on the light cone. For the causal line integrals, we use 
the short notation 




A.l Electromagnetic Potential 

AP{x,y) = jy^eiz-- 

iiOTT-' Jx 

+ - 6a^ + 2a) ^ 7' (nF^j) logz 
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(\a^-a) e'^^' {uFi^)^klH log- 

Jx 
ry 
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Jx 



m / Fij cr*"' log z 

Jx 

32.3 ■.J\a'-a)j,eios. 



647r3 
ie 
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m 
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647r3 

ie 

~ 1287r3 



m 



ry 

Jx 

r{2a-l)j'Fije log; 

' r e'^"' F,,^,^^^i log, 

Jx 



y 
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with the electromagnetic field tensor Fjk = djA^ — d^Aj and the electromagnetic current 
A. 2 Axial Potential 
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with the axial field tensor Fjk = djE^ — dkBj and the axial current = diF^K 

A. 3 Scalar Potential 
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A. 4 Pseudoscalar Potential 
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A. 5 Bilinecir Potential 
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B Perturbation by a Gravitational Field 

In this appendix, we outline how the light-cone expansion can be extended to a perturbation 
by a gravitational field. For the metric, we consider a perturbation hjk of the Minkowski 
metric rfj^ = diag(l, —1, —1, —1), 

gjk{x) = rjjk + hjk{x) . 

We describe gravitation with the linearized Einstein equations (see e.g. [6]). According to 
the usual formalism, we raise and lower tensor indices with respect to the Minkowski met- 
ric. Using the transformation of hjk under infinitesimal coordinate transformations, we can 
assume Par. 105] that 



d'^hjk = - djh with h := h'^ 



In the so-called symmetric gauge, the Dirac operator takes the form [7] 

In contrast to ()1.2p . the perturbation is now itself a differential operator. 

One complication arises from the fact that the integration measure in curved space is 
yj\g\ d^x = (\ + d^x, whereas the formula (jl.7p for the perturbation of the Dirac sea is 



2' 

valid only if one has the integration measure d'^x of Minkowski space. Therefore we first 
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transform the system such that the integration measure becomes d^x, then apply p.7p . and 
finally transform back to the original integration measure ^/\g\ d^x. Since the scalar product 

y^$^(i^x = J ilgl'i'^) {\g\^ ^) d^x 

is coordinate- invariant, the transformation to the measure d^x is accomplished by 

The perturbation of the transformed system is given by ()1.7p . 

AP('^'-)(x,y) = -jd^z(^s{x,z)(^-'-^^h]^-'-{m^))P{^.y) 

+ P{x,z)^-'-j^h^^-'-{mz))s{z,y)^ . (B.l) 

The formula for the transformation of the Dirac sea to the original integration measure 
Vldl d'^x is 

Pix,y) + AP{x,y) = \g\-Hx) \g\-"^ (y) [Pix,y) + Ap(^'")(x, y)) . 

Thus 

APix,y) = Ap('^'-)(x,y) - ^ + %)) P(x, y) . 

The factors P{z,y) and s{z,y) in (jB.ip depend only on (z — y), i.e. 

so we may rewrite the z-derivatives as y-derivatives, which can be pulled out of the integral. 
Furthermore, the relations 

J d^z P{x, z) {i^,h{z)) s{z, y) = J d^z P{x, z) - m), h{z)\) s[z, y) 

= -P{x,y)h{y) 
J d^z s{x, z) {i^M^)) P{z, y) = h{x) Pix, y) 

allow us to simplify the factors {^h) in the integral. In the resulting formula for AP{x,y), 
one recovers the perturbation by an electromagnetic potential. More precisely, 

AP(x,y) = (-Ihix) - ^%))p(x,y) - ^AP[^^ h';]ix,y) , (B.2) 

where APlj^hjjlx, y) is the perturbation (jl.7p of the Dirac sea corresponding to the electro- 
magnetic potential B = j^hj. The light-cone expansion of AP{x, y) is obtained by substitut- 
ing the light-cone expansion of AP[7-'7i^](x, y) into (|B.2p and calculating the y-derivatives. 
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To the order 0{{y — x)^) on the hght cone, this gives the following formula for the light-cone 
expansion of the Dirac sea in the gravitational field: 
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where Rjk and are the (linearized) Ricci tensor and scalar curvature, respectively. A general 
difference to the formulas of Appendix |A] is that AP{x,y) now has a stronger singularity on 
the light cone. This is a consequence of the y-derivative in ()B.2p . The leading singularity 
of AP{x, y) can be understood as describing the "deformation" of the light cone by the 
gravitational field in linear approximation. 

We finally remark that this method works also for the higher-order perturbation theory as 
developed in [5]. It can be used to perform the light-cone expansion of higher-order Feynman 
diagrams in the presence of a gravitational field. 
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